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Abstract 

MacMahon proved a simple product formula for the generating function of the plane 
partitions fitting in a given rectangular box. The theorem implies the number of lozenge 
tilings of a semi-regular hexagon on the triangular lattice. By investigating the lozenge 
tilings of a hexagon with a hole on the boundary, we generalize the ordinary plane partitions 
to piles of unit cubes fitting in a union of several adjacent rectangular boxes. We extend 
MacMahon’s classical theorem by proving that the generating function of the generalized 
plane partitions is given by a simple product formula. 

Keywords: perfect matchings, plane partitions, lozenge tilings, dual graph, graphical con¬ 
densation. 


1 Introduction and main results 


Given k positive integers Ai ^ A2 ^ ^ A A plane partition of shape (Ai, A2, • • •, A j~) is an 

array of non-negative integers 


^ 1,1 n p 2 n 1,3 . ... 

^2,1 ^2,2 ^2,3 . ™2,A 2 


nk, 1 nk, 2 nk,3 ■ ■ ■ n k,x k 

so that n l j ^ fbj+i and riij ^ (i.e. all rows and all columns are weakly decreasing from 

left to right and from top to bottom, respectively). The sum of all entries of a plane partition 
7T is called the volume (or the norm ) of the plane partition, and denoted by |7r|. 

The plane partitions of rectangular shape (b, b, ..., b ) (a rows) with entries at most c are 
usually identified with their 3-D interpretations — piles (or stacks) of unit cubes fitting in an 
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by the National Science Foundation. 


1 



a=2 


Figure 1.1: The shamrock 54 , 2 , 2 , 3 - 


a x b x c box. The latter are in bijection with lozenge tilings of a semi-regular hexagon of side- 
lengths a,b,c,a,b,c (in clockwise order; starting from the northwest side 1 ) on the triangular 
lattice. Here, a lozenge (or unit rhombus ) is union of any two unit equilateral triangles sharing 
an edge; and a lozenge tiling of a region is a covering of the region by lozenges so that there are 
no gaps or overlaps. 

Let q be an indeterminate. The q-integer [n] q is defined by [n] q := 1 + q + q 2 + ... + q n ~ l . 
We also define the q-factorial by [n] q \ := [l] g [2] g ... [a] 9 , and the q-hyperfactorial function by 
H g (n) := [0]^!- [l]g!• [2] g !... [n—1] 9 !. MacMahon [Ma] proved that the volume generating function 
of the plane partitions fitting in an a x b x c box is given by 

| 7 r| _ Hq(a) Hq,(ft) H g (c) H g (q + b + c) 

2^ Hq(q + b) Hq(6 + c) H g (c + a) 


The q = 1 specialization of the Macmahon formula (1.1) is equivalent to the fact that the number 
of lozenge tilings of a semi-regular hexagon Hex(a, b, c ) of side-lengths a, b, c, a, b, c is equal to 

H(q)H(fr)H(c)H(q + fr + c) 

H(q + b) H(6 + c) H(c + q) ’ { ' ) 

where H(n) := Hi(n) = 0! • 1! • 2!... (n — 1)! is the ordinary hyperfactorial function. 

Extending the MacMahon’s classical theorem (for the case when q = 1), Ciucu, Eisenkolbl, 
Krattenthaler, and Zare [CEKZ] proved a simple product formula for the number of tilings of 
a hexagon of side-lengths a,b + m,c,a + m,b,c + m with a triangular hole of size m in the 
center (the region was called cored hexagon in [CK]). Recently, Ciucu and Krattenthaler [CK] 
generalized the latter result further by extending the central triangular hole to a shamrock hole 
consisting of four adjacent equilateral triangles, and showing that the tiling formula is still a 
simple product formula. Precisely, the shamrock S rn a i- /c is the union of four equilateral triangles 
with sides m, q, b, c on the triangular lattice described in Figure 1.1. 

Let us consider a related situation of the cored hexagon in [CEKZ] when the triangular hole 
appears on the boundary (instead of in the center). It has been proven that the regions of such 
type have the number of lozenge tilings given by a simple product formula (see e.g. Proposition 
2.1 in [CLP]). As suggested by Ciucu and Krattenthaler’s work in [CK], we extend the triangular 
hole on the boundary of the hexagon to a shamrock hole as follows. 

1 From now on, we always list the side-lengths of a hexagon on the triangular lattice in the clockwise order, 

starting from the northwest side. 
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x+c y+b 

Figure 1.2: Hexagon with a shamrock removed along the boundary. 


We start with a hexagon of side-lengths z + a + b + c, x + y + m, t + a + b + c, z + m,x + 
y + a + b + c,t + m. Next, we remove a shamrock S ma! i, !C from the base of the hexagon so that 
the lower-left vertex of the a-triangle in the shamrock is x + c units to the right of the lower-left 

I X V z t\ 

vertex of the hexagon. We denote by Q[ , the resulting region. Figure 1.2 shows 

\m a b c) 

the region Q 

Interestingly, the lozenge tilings of our new region are always enumerated by a simple product 
formula. 


(4 3 2 3\ 
\4 3 2 1 )■ 


Theorem 1.1. For non-negative integers x, y, z , t, m, a, b, c 


M (o( X y Z ^ ^ — H (m + a + b + c + x + y + z + t) 

V V m a b C J J H(m + a + b + c + x + y + + a + b + c + x + y + z) 

H (m + a + b + c + x + t) H (m + a + b + c + x + y) H(m + a + b + c + y + z) 
H(m + a + b + c + z + t ) H(m + a + b + c + x ) H(m + a + b + c + y) 

H(a) H (y) H(z) H(t) H(m) 3 H(q) 2 H(b) H(c) H (m + a + b + c) 

H(x + t) H(y + z) H(m + a) 2 H(m + b) H (m + c) 

H(m + b + c+z + t ) H(m + a + c + x) H [m + a + b + y) 

H(m + b + y + z) H (m + c + x + t) 

H(c + x + t) H(6 + y + z) 

H(a + c + x) H(a + b + y) H (b + c + z + t) ’ 

where we use the notation M (R) for the number of lozenge tilings of a region R 


Next, we consider a g-analog of Theorem 1.1. 

Similar to the bijection between lozenge tilings of a semi-regular hexagon Hex(a,b, c) and 

f x y z 

plane partitions fitting in an a x b x c box, one can view a lozenge tiling of Q ( ' ^ ^ J as 
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(a+b+c-m,t+m) 


3 

(a+b,t+c) 


_ 


2 

(0,1+a+b+c) 


1 

<0.a) 


y+b 


(C) 


Figure 1.3: (a) Viewing a lozenge tiling of the region as a pile of unit cubes in a special box. 
(b) The tiling corresponding to the empty pile, (c) The projection of the compound box B on 
the Oij plane. 


( X V z t\ 

7 , which is the union of 6 

m a b c) 

adjacent rectangular boxes (see Figure 1.3(a) for x = y = 2 = 3, t = 4, m = 3, a = b = c = 2). 
Figure 1.3(b) gives a 3-D picture of the compound box B by showing the empty pile; the bases 
of the rectangular boxes in B are labeled by 1, 2 ,..., 6. 

Projecting the compound box B on the Oij plane, we get a projective diagram as in Figure 
1.3(c). In this diagram, each rectangular box in B is represented by a rectangle with a pair of 
integers (a, 6), where a is the level of the base and b is the height of the box. We always assume 
that the base of the box 1 is on level 0. We also note that the rectangles corresponding to the 
boxes 1 and 3 are overlapped (the intersection is indicated by the shaded area in Figure 1.3(c)). 
However, these two boxes are not overlapped since the base of the box 3 is above the top of the 
box 1 (as a + b ^ a). 

We call the piles of unit cubes fitting in the compound box B generalized plane partitions , 
since they have a similar monotonicity as (the 3D-interpretation of) the ordinary plane partitions: 
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the tops of their columns (of unit cubes) are weakly decreasing along Oi and Oj. 

Similar to MacMahon’s theorem (1.1), we have a closed form product formula for the volume 
generating function of the generalized plane partitions. 


Theorem 1.2. Let m, a, b, c, x, y, z, t be non-negative integers. Then 


X> H 


Hg (m + a + b + c + x + y + z + t) 

Hg (m + a + b + c + x + y + t ) Hg(m + a + b + c + x + y + z ) 

Hg(m + a + b + c + x + t) H q [rn + a + b + c + x + y) Hg(m + a + b + c + y + z ) 

Hg(m -\-a + b + c + z + t ) H q {m + a + b + c + x) Hg (m + a + b + c + y) 

Hg(x) Hg ( 'y ) Hg(z) Hg (t) Hg(m ) 3 Hg(a ) 2 Hg(b) Hg(c) Hg(m + a + b + c) 

Hg(x + f) Hg(y + z) Hg (m + a) 2 Hg(m + b ) Hg(m + c) 

Hg(m + b + c + 0 + t) Hg(m + a + c + x) H g (m + a + b + y) 

Hg(m + 6 + y + z) Hg(m + c + x + t) 

_Hg(c + x + t) Hg (b + y + z )_ 

Hg(a + c + x) Hg(a + b + y) H q (b + c + z + t) ’ 


where the sum, on the left-hand side is taken over all generalized plane partitions n fitting in the 
x y z fi 


compound box B 
in it). 


m a 


, and where |7r| is the volume of n (i.e. the number of unit cubes 


One readily sees that Theorem 1.1 is a consequence of Theorem 1.2. Moreover, by letting 
m = a = b = c = 0 (i.e. the shamrock hole is empty), we get the MacMahon’s classical theorem 
(1.1) from Theorem 1.2. 

The goal of the present paper is to prove Theorem 1.2 by using the graphical condensation 
method first introduced by Eric H. Kuo in [Ku04j. This condensation can be viewed as a combi¬ 
natorial interpretation of Dodgson condensation (also called as the Jacobi-Desnannot identity, 
see e.g. [Mu], pp. 136-149). We refer the reader to e.g. [YYZ], [YZ], [Ku06], [Sp], [Ci], [Fu] for 
various aspects and generalizations of the method; and e.g. [CK], [CL], [CF14], [CF15], [KW], 
[Lal5a], [Lal5b], [Lal5c], [LMNT], [Zh] for recent applications of Kuo condensation. 

The rest of our paper is organized as follows. For ease of reference, we quote the several 
preliminary results in Section 2, including the particular version of Kuo condensation employed 
in our proofs. In order to apply Kuo condensation to our Q-type regions, we consider several 
simple weight assignments on the lozenges of the regions in Section 3. Section 4 generalizes a 
related work of Ciucu and Krattenthaler in [CK, Theorem 3.1] by (/-enumerating lozenge tilings 
of a magnet bar region (the b = c = 0 specialization of a Q-type region). Finally, we prove 
Theorem 1.2 in Section 5. 


2 Preliminaries 

Let G be a finite simple graph without loops. A perfect matching of G is a collection of disjoint 
edges covering all vertices of G. Let R be a region 2 . The (planar) dual graph of R is the graph 
whose vertices are unit triangles in R and whose edges connect precisely two unit triangles 

2 From now on, we use the word region to mean a finite connected region on the triangular lattice. 
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sharing an edge. One can identify the lozenge tilings of R with the perfect matchings of its dual 
graph. 

For a weighted graph G, we define the matching generating function M(G) of G to be the 
sum of the weights of all perfect matchings in G, where the weight of a perfect matching is the 
product of weights of its edges. If the lozenges of a region R are weighted, we define similarly 
the tiling generating function M (R) of R. We notice that each edge of the dual graph G of the 
region R carries the same weight as its corresponding lozenge in R. 

The following Kuo’s condensation theorem is the key of our proofs. 

Theorem 2.1 (Theorem 5.1 in [Ku04]). Let G = (Vi,V 2 ,E) be a (weighted) bipartite planar 
graph in which |Vi| = \Vf\- Assume that u,v,w,s are four vertices appearing in a cyclic order 
on a face of G so that u,w £V± and v,s € Vi- Then 

M(G) M(G — {u, v. w, s}) = M(G — {u,«}) M(G — {re, s}) + M(G — {u, s}) M(G — {v, w}). (2.1) 

A forced lozenge of a region I? is a lozenge contained in any tilings of R. Assume that we 
removed several forced lozenges £ 1 ,( 2 , ■ ■ ■ An from the region R, and denote by R' the resulting 
region. Then one clearly has 

n 

M(i2) = M(i2 , )IJ^i), (2.2) 

i= 1 

where wt(£i ) denotes the weight of the lozenge l{. 

If a region R admits a lozenge tiling, then the numbers of up-pointing unit triangles and 
down-pointing unit triangles in R are equal. Moreover, if a region satisfies the latter balancing 
condition, we say that the region is balanced. The following lemma can be considered as a 
generalization of the identity (2.2). 

Lemma 2.2 (Region-splitting Lemma). Let R be a balanced region. Assume that a subregion S 
of R satisfies following two conditions: 

(i) (Separating Condition) There are only one type of unit triangles (up-pointing or down¬ 
pointing) running along each side of the border between S and R — S. 

(ii) (Balancing Condition) S is balanced. 

Then 

M(R) = M{S)M(R- S). (2.3) 

Proof. Let G be the dual graph of the region R. Then the dual graph K of the subregion S is 
an induced subgraph of G. It is easy to see that K satisfies the conditions in Lemma 3.6(a) in 
[Lal4], and the lemma follows. □ 

3 g-weight assignments 

Lozenges in a region R come with three different orientations: left, right, and vertical lozenges 
(see Figure 3.1). Next, we consider three simple g-weight assignments of lozenges in our region 

Q:=q( X y Z ^ as follows. 

\m a b cl 
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Vertical 


Right 


Left 


Figure 3.1: Three orientations of lozenges. 


View a lozenge tiling T of the region Q as a generalized plane partition (i.e. a pile of unit 
cubes); each right lozenge is viewed as the top of a column of unit cubes. We assign to each right 
lozenge a weight q x , where x is the number of unit cubes in the corresponding column. All left 
and vertical lozenge are weighted by 1. We call this assignment the natural q-weight assignment 
of Q , and use the notation wto for the assignment (see Figure 3.2(a) for x = y = 2 = 3, t = 4, 
m = 3, a = b = c = 2). 

Besides the natural q- weight assignment wto, we consider the following two simple weight 
assignments: 

(1) Assignment 1. The weights of left and vertical lozenges are all 1. The weight of a right 
lozenge is q l , where l is the distance between the left side of the lozenge and the southeast 
side of the region Q. We use notation wt\ for this weight assignment (see Figure 3.2(b)). 

(2) Assignment 2. All left and vertical lozenges are still weighted by 1. However, a right 
lozenge is now weighted by q n , where n is the distance between the top of the lozenge and 
the base of the region Q. This assignment is denoted by wt2 (see Figure 3.2(c)). 

Let T be a tiling of Q. We denote by wto(T), wt\(T ) and wt2{T) the weights of the tiling 
T with respect to the weight assignments wto, wt\ and wt2- We also denote by Mo(Q), Mi(Q) 
and M2(Q) the tiling generating functions of Q corresponding to the weight assignments wto, 
wt\ and wt2- It is easy to see that Mo(Q) is exactly the volume generation function of the 
generalized plane partitions corresponding to the lozenge tilings of Q. 

In the next proposition, we will show that the three weight assignments are the same up to 
some multiplicative factors. 

We define two functions 


x y z t 
m a b c 


:= m 


y + b+l 
2 




y +1 


+ m(z + b) (y + a + b) + (z + b) 


+ x(z + b + c)(y + m + a + b + c) + (z + b + c) 


x + 1 
2 


+ a{x + c)(y + a + b) + a 


m + 1 
2 

x + c + 1 
2 


and 


g 


x y z t 

m a b c 

z + b+l 
2 


:= (y + b) 


m + 1 
2 


+ myz + y 


+ m 


+ x(m + a)(z + b + c) + x 


z + 1 
2 

z + b + c + 1 
2 


+ m(z + b)(m + a) 
+ (x + c) 


a + 1 
2 


(3.1) 


(3.2) 
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Figure 3.2: Three g-weight assignments on a sample tiling of region Q: (a) wto, (b) wt±, (c) wt 2 - 
The right lozenges with label x are weighted by q x . 
























































































Proposition 3.1. For any non-negative integers m,a,b,c,x,y, z,t 


Mi Q 


and 


M 2 Q 


x y z t 
m a b c 


x y z t 
m a b c 


= q 


= q 


x y z t 
m a b c 


x y z t 
m a b c 


T.F 


TF. 


(3.3) 


(3.4) 


where the sums on the right-hand sides are taken over all generalized plane partitions it fitting 

x y z tfi 


in the compound box B 


m a b c 


Proof. We use the following shorthand notations in this proof: f := f 

%{* y l T B-.= B ( X « l ‘land Q :=Q 

\m a b cl \m a b cl 


g : = 


x y z t 
m a b cj 

x y z t\ 
m a b cj ' 

Let T be any lozenge tiling of the region Q, and n the generalized plane partition corre¬ 
sponding to T. We only need to show that 


wti (T) _ f 


wt 2 (T) 


q 


M 


and —= (7 s . 

aM 


(3.5) 


Assume that the box i of the compound box B has size a* x bi x c, (for 1 ^ i ^ 6 ). The base 
of the box i is pictured as a parallelogram T) in Figure 1.3(b). We assume in addition that the 
left side of Pi is Xj units to the left of the southeast side of the region Q, and the bottom of Pi 
is yi units above the bottom of region Q. 

Divide the generalized plane partition 7 r into 6 disjoint sub-partitions 7 Tj (1 ^ i ^ 6 ) fitting in 
the box i. Each partial-partition vr^ in turn gives a lozenge tiling Tj of the semi-regular hexagon 
Hex(ai,bi,Ci). Figure 3.3(a) shows the partial-partition 713 of the generalized plane partition 7 r 
in Figure 1.3(a), as well as, the relative positions of the parallelogram P 3 to the bottom and the 
southeast side of the region Q (for x = 3, y = 3, z = 3, t = 4, m = 3, a = 2, b = 2, c = 2). 

Apply the weight assignment wt\ to the whole tiling T of the region Q. This yields a local 
weight assignment wt^ for the tiling T t of hexagon Hex(ai,bi, cfi). Precisely, each right lozenge 
in Tj is now weighted by q Xi+l l where l is the distance between the left side of the lozenge and 
the southeast side of the hexagon Hex(ai,bi,Ci). Encode the tiling Tj as a family of bi disjoint 
lozenge-paths connecting the top and the bottom of the hexagon (see the dotted paths in Figure 
3.3(b)). Dividing the weight of each right lozenge in the lozenge-path j (from right to left) by 
q x i+fi we get the weight assignment wto for Tj. Since the end points of the above lozenge-paths 
are fixed, each lozenge-path has exactly Oj right lozenges. Thus, we have 


wto{Ti ) 


fll 7r i| 


_ ^aibiXi+aibi{bi+ 1)/2 


Multiplying all above equations for i = 1,2,..., 6, we get 


wti(T) 


— qJ2^=l a i b i x i+ a i b i(. b i + 1 )/ 2 


(3.6) 
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(a) 



Figure 3.3: The partial-partition corresponding to the room 3. 
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Next, we assume that the whole tiling T of Q is weighted by wt 2 - We now encode the tiling T* 
of Hex(cii,bi, Ci) as an a*-tuple of disjoint lozenge-paths connecting the northwest side and the 
southeast side of the hexagon (see Figure 3.3(c)). Dividing each right lozenge in the lozenge-path 
j (from bottom to top) by q yi+J , we get back again the weight assignment wto. We note that 
each lozenge-path has now bi right lozenges. Similar to the case of wt\, we have 


wt '2 (T) _ J2i=i aibiyi+biCLifai+l) / 2 


(3.7) 


Obtaining the formulas of Oj, bi, Xi, y % in terms of m, a, b, c, x, y, z, t from Figures 1.3(b) and (c), 
we get f = Ya=i a ibiXi + a,ibi{bi + l)/2 and g = )T^ =1 + ^(a, + l)/2. This finishes our 

proof. □ 


We note that the powers q f and q s in the above proposition are exactly the weights wti(To) 
and wt 2 (To) of the tiling To corresponding to the empty pile of unit cubes in Figure 1.3(b). 

x y z t'' 
y m a b Cj 

empty shamrock hole. We have the following consequence of Proposition 3.1 and MacMahon’s 
q- formula (1.1). 


View the hexagon Hex(a, b, c ) as a special case of the region Q := Q 


with an 


Corollary 3.2. For non-negative integers a,b,c 


and 


Mi ( Hex(a, b, 


\\ _ ab(b+ 1)/2 Hg( a ) Hg(fr) H g (c) H g (a + b + c) 
H 9 (a + b) H g (6 + c) H ? (c + a) 


M 2 ( Hex(a,b,c )) = q ba ^+ 1 )/ 2 


H,(a) H q (b) H g (c) H g (a + 6 + c) 
H q (a + b) H g (6 + c) H g (c + a) 


(3.8) 


(3.9) 


A column-strict plane partition is a plane partition having columns strictly decreasing. A 
semihexagon SH a i, is the upper half of a lozenge hexagon Hex(a, b, a). We are interested in the 
lozenge tilings of the semihexagon SH a ^, where a up-pointing unit triangles at the positions 
si,S 2 , ■ ■ ■ ,s a have been removed from the base. Denote by SH a j,(si, s 2) • • •, s a ) the resulting 
semihexagon with dents (see Figure 3.4(a) for the region 2, 6, 7,10,11,12)). Assume 

that the lozenges in the semihexagon are weighted by wt 2 , and we still use the notation M 2 
for the corresponding tiling generating function of the semihexagon with dents. There is a 
well-known (weight preserving) bijection between the lozenges of SH a ^(si, s 2 ,..., s a ) and the 
column-strict plane partitions of shape (s a — a, s a _i — a + 1,... , si — 1) with positive entries at 
most a, i.e. ief 2 (T) = q^ nT ^, where ttt is the plane partition corresponding to the tiling T (see 
e.g. [CLP] and [CS]). 

We have the following ^-enumeration of the lozenge tilings of a hexagon with a triangular 
hole on the base K a (x,y,z,t ) (defined as the region restricted by the bold contour in Figure 
3.4(b)). 


Lemma 3.3. For non-negative a, x, y, z, t 


™ (ts .. _ ^ H q {a)R q (x)Jl q (y)Jl q (z)R q (t) 

M 2 (K a (x,y,z,t)) -q y y 2 > 1 TTTTvfr T ^ rr u. 


Hg(x + t) H q (a + x) H q (a + y) H q (y + z) 

H g (a + x + t) H g (a + x + y) H q (a + y + z) H 9 (a + x-\-y + z-\-t) 
H q (a + x + y + t) H q (a + x + y + z) H 9 (a + t + z) 


(3.10) 
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b=5 


x+y 




(b) 


Figure 3.4: (a) The semihexagon with dents SH? 1, 2, 6, 7,10,11,12). (b) Obtaining the region 
K a (x,y,z,t ) (restricted by the bold contour) from a semihexagon with dents. 


x+y+m 



Figure 4.1: The magnet bar i? 2 , 2 (4,3, 3, 2). 

Proof. By the above bijection between lozenge tilings of the semihexagon and the column-strict 
plane partitions, we have 


M 2 (SH aib ( ai ,S 2 , • • •, s a )) = J2 = q ES = l(81 "° 






q j - q l ’ 


(3.11) 


where the sum after the first equality sign is taken over all column-strict plane partitions y of 
shape (s a — a, s a ~ i — a + 1,..., si — 1) with positive entries at most a, for the second equality 
see e.g. [St, pp. 374-375]. 

The region K a (x,y, z,t) is obtained by removing forced vertical lozenges from the serni- 
hexagon SH a+z+ttX + y with dents at the positions {1,2, ,.,t} U {t + x + l,t + x + 2,... ,t + x + 
a} \J {t + x + a + y + l,t + x + a + y + 2,... ,t + x + a + y + z}. Thus, our lemma follows from 

(3.11). □ 


4 Two g-enumerations of magnet bar regions 

I X V z t\ 

When b = c = 0, our region Q , becomes a magnet bar region B ma (x,y, z,t) first 

a b cj 

introduced in [CK], Figure 4.1 shows the magnet bar region i?2,2(4, 3, 3,2). In this section, we 
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(q-) enumerate lozenge tilings of the magnet bar region B m ^ a {x,y,z,t). Our proof of Theorem 
1.2 will use these ^-enumerations. 

Proposition 4 . 1 . For non-negative integers m,a,x,y, z,t 

m 2 (B mia (x,y,z,t)) = q y( rn 2 1 )+( m + x +y)( z t 1 )+ m y z +( m + a )( x + m ) z M a V) 

H q (m +a+x+y+z+t ) 

H g (m + a + x + y + t) H q (m + a + x + y + z ) 

H g (m + a + x + t) H q (m + a + x + y) H q (m + a + y + z) 

H q (m + a + z + t) H q (m + a + x) H q (m + a + y) 
y H g (x) H q (y) H q (z) H q (t) H q (m) H q (a ) 2 

H g (a + x) H g (a + y) H q (z + t) H q {m + a) 

R q (m + z + t) R q (m + a + x) R q (m + a + y) 

H q (m + y + z) H q (m + x + t) 

Proof. We prove (4.1) by induction on y + z + 1. Our base cases are the situations when m = 0, 
a = 0, y = 0, z = 0 or t = 0. 

Throughout this proof, we assume that our magnet bar region is weighted by wt2- 
If m = 0, then our magnet bar region becomes the region K a (x,y, z,t) in Lemma 3.3, and 
(4.1) follows. 

If a = 0, by removing forced lozenges along the base of the region B m fi(x,y, z,t), we get 
the weighted hexagon Hex(z , x + y + m, t) in which a right lozenge is weighted by q m+l , where 
l is the distance from the top of the lozenge to the bottom of the hexagon (see Figure 4.2(e)). 
By dividing the weight of each right lozenge of the hexagon by q m , we get back the weight 
assignment wt2- Since the product of weights of the forced lozenges in Figure 4.2(e) is equal to 

/ m+l\ 

qy\ 2 i )We get from (2.2) 

M 2 (fl m , o(x, y, z, t)) = qy^q^+y+m) M2 (^ ex(Zj x + y + m , t j), (4.2) 

where the factor q mz ( x +y+ m ) comes from the weight division. Then (4.1) follows from Corollary 
3.2. 

If y = 0, after removing forced vertical lozenges (which have the weight 1), we get a new 
weighted region R (the region restricted by the bold contour in Figure 4.2(b)). By rotating R 
60° clockwise and reflecting it about a vertical line, we get the region K m (z,a,x,t ) weighted by 
wt±. Thus, we have 

M 2 (B mt a(x, 0, 2 , t)) = Ml ( K m {z , a, x, t)) , (4.3) 

and (4.1) follows from Proposition 3.1 and Lemma 3.3 (one can view any LT-type region as the 
m = b = c = 0 specialization of a Q-type region). 

If z = 0, by applying Region-splitting Lemma 2.2 as in Figure 4.2(c), we get 

M 2 (B mta (x,y, 0 ,t)) = M 2 (Hex(m,y,a))M 2 ^B m ^ a (x,y, 0 ,t) - Hex(m,y,a)j. (4.4) 

Next, we remove the forced left lozenges from the region B m>a (x, y, 0, t) — Hex(m,y,a ) and 
obtain the hexagon Hex(a,x,t + m) weighted by rct 2 . Then we get 

M 2 {B mta (x,y, 0 ,t)) =M 2 {Hex(m,y,a))M2(Hex(a,x,t + m)), (4.5) 
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y+m 


x+m 




* y +m x+v+m 



Figure 4.2: The base cases in the proofs of Propositions 4.1 and 4.2: (a) x = 0, (b) y = 0, (c) 
z = 0, (d) t = 0, and (e) a = 0. 


and (4.1) follows from Corollary 3.2. 

If t = 0, similar to the case when z = 0, the Region-splitting Lemma 2.2 implies 

M 2 y, z , 0)) = M 2 ( Hex[z + m, y, a)) M 2 (^B m , a [x, y, z, 0) - Hex[z + m, y, a)) (4.6) 

(see Figure 4.2(d)). We also get the hexagon Hex[a,x,m ) (weighted by iut 2 ) after removing 
forced lozenges from the region B m ^ a {x, y, z, 0) — Hex(z + m, y, a). However, our forced lozenges 
are now right lozenges, which have weight product equal to g( m + a )( 3: + m ) 2 +( :r + m )( 2 ). Thus, we 
get 


M 2 (^B m ^ a (x, y, z, 0) - Hex(z + m, y, a)) = g,( m + a )( 3: + m ) 2 +( a: + m )( t ) m 2 [Hex(a, x, m )), 
so 

M 2 ( B mta (x , y, z, 0)) = q( m + a )( x + m ) z +( x + m ){ t ) m 2 ( Hex(z + m, y, a)) M 2 [Hex(a, x, m)). 

(4.7) 

Again, (4.1) is implied by Corollary 3.2. 

For the induction step, we assume that m,a,y,z,t ^ 1 and that (4.1) holds for any magnet 
bar regions, which have the sum of the y-, z- and t-parameters strictly less than y + z + t. 

We apply Kuo Theorem 2.1 to the dual graph G of the magnet bar region B miCL [x,y,z,t) 
(weighted by rcf 2 ). We pick the four vertices u, v, w, s as in Figure 4.3(b). In particular, the four 
shaded unit triangles correspond to the four vertices: the shaded unit triangle corresponding 
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x+y+m 


x+y+m 



Figure 4.3: Obtaining the recurrence for the tiling generating functions of magnet bar regions. 
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to u is the lowest one, and v, w, s correspond to the next shaded unit triangles as we move 
counter-clockwise from the lowest one. We notice that the north side of the region has length 
x + y + m^-y + m^-2 and the northeast side has length t + a ^ 2, so the four vertices u, v, w, s 
are well-defined. 

By removing the lozenges forced by the shaded unit triangles, we get back new £>-type regions 
weighted by ref 2 - Collecting the weights of those forced lozenges, we get 

M(G - {U, V, W, s}) = q (* + ™ +1 )H*+y+rn-2)( Z+ t + m + a) ^ (^(^ y _ W - 1)) , (4.8) 


/z+ra+l\ , v 

M(G - {«,u}) = 2 > M 2 {B m ^ a (x, y - 1, 2 , t)) , (4.9) 

M(G - {w, s}) = q ^+y+m-2)(z+t+m+a) M2 y,z,t- 1)) , (4.10) 

M (G-{u,s}) = q( 2 )M 2 (B mt a(x,y -l,z + l,t-l)), (4.11) 

and 

M(G - {u, re}) = q ^+y+m-i)(z+t+m+a) ( Sm a ( Xj y :Z _ i > ()) (4.12) 


(see Figures 4.3(b)-(f), respectively). Plugging the above identities into the equation (2.1) in 
Kuo Condensation Theorem 2.1, we obtain 

M 2 (B mia (x, y, z, t)) M 2 (B m>a (x, y-l,z,t-l)) = 

M 2 (B m ^ a (x, y-l,z, t)) M 2 (B miCL (x, y,z,t- 1)) 

+ q z+t+m+a M 2 (B m , a (x, y - 1, z + 1, f - 1)) M 2 ( B m} a(x , y,z-l,t)). 

(4.13) 

All regions in the above equation, except for the first one, have the sum of their y-, z- and 
f-parameters strictly less than y + z + f. Thus, by the induction hypothesis, those regions have 
their tiling generating functions given by (4.1). By substituting these formulas into the above 
equation and working on some simplifications, one readily gets M 2 (B mtCL (x, y, z, t)) equal exactly 
to the expression on the right-hand side of (4.1). This finishes our proof. □ 

We need another (/-enumeration of lozenge tilings of the magnet bar region as follows. 

Assume that we now give all right and left lozenges in the magnet bar region B m ^ a (x, y,z,t) 
a weight 1. Next, we give a vertical lozenge a weight q l , where l is the distance between the 
northeast side of the lozenge and the southwest side of the region. We denote by wt% the new 
weight assignment, and M 3 the corresponding tiling generating function. 

Proposition 4.2. For non-negative integers m,a,x,y, z,t 

M 3 {B m j x , y ,z,t)) = ,™("r)+‘(‘ + r‘)+«(*+™><-+“>+«(‘ + r‘) 

H g(m + a + x + y + z + t) 

H q (m + a + x + y + t ) H q (m + a + x + y + z) 

H q (m + a + x + t) H q (m + a + x + y) H q (m + a + y + z) 

B q (m + a + z + t) H q (m + a + x) H q (m + a + y) 
x HgQc) H q (y) H q (z) H q (t) H q (m) H q (a ) 2 

H g (a + x) H q (a + y) H q (z + f) H q (m + a) 

H q (m + z + t) H q (m + 0 + x) H q (m + a + y) (a-\a'\ 

H q (m + y + z) H q (m + x + t) 
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Proof. The equality (4.14) can be treated similarly to (4.1) in the Proposition 4.1 by induction 
on y + z + t. The base cases are the situations when a = 0, x = 0, y = 0, z = 0 or t = 0. 

Assume that our magnet bar region is weighted by wt^ 3 . 

If a = 0, we get a weighted version of the hexagon Hex(z : x + y + m,t) by removing forced 
lozenges from the magnet bar reigon as in Figure 4.2(e). Rotating the hexagon 60° clockwise 
and reflecting it over a vertical line, we get the hexagon Hex(z,t,x + y + m) weighted by wt 2 , 
and (4.14) follows from Corolarry 3.2. 

If x = 0, after removing forced vertical lozenges (whose weight product is equal to 2 )) 

as in Figure 4.2(a), we get a weighted region R. Next, we rotate R 60° counter-clockwise and 
reflect it about a vertical line. This way, we get the weighted region K m (a,t, z,y) in which a 
right lozenge is weighted by q a+l , where l is the distance from the top of the lozenge to the 
bottom of the region. We divide the weight of each right lozenge in the latter region by q a , and 
get back the weight assignment wt 2 - Thus, 

M 3 {B mia (0,y,Z,t)) = q (t+m)(^) q a Z t+aHz + m) M2 (4.15) 

where the factor q azt + a2 ( z + m ) comes from the weight division. Then (4.14) follows from Lemma 
3.3. 

If y = 0, by removing forced lozenges (whose weight product is q( x + a )( z + m ) a + a ( 2 )) and 

rotating the resulting region 60° clockwise, we get a weighted version of region K m (a, z,t,x) 
in which a right lozenge is weighted by q m + a + x + z + 1 ~f where l is the distance from the left 
side of the lozenge and the southeast side of the region (see Figure 4.2(b)). By dividing the 
weight of each right lozenge by q m + a + x + z + l ; we get back the weight assignment wt\, where q 
is replaced by <7 . Thus, (4.14) follows from Proposition 3.1, Lemma 3.3 and the simple fact 

[n\ q -i = [nlq/q^-^. 

If z = 0, we apply Region-splitting Lemma 2.2 (and remove forced lozenges weighted by 1) 
to split our region into two hexagons as in Figure 4.2(c). Next, we rotate the right hexagon 
60° counter-clockwise and reflect it about a vertical line to get the hexagon Hex(a,t + m,x) 
weighted by wt 2 - For the left hexagon, we also rotate it 60° counter-clockwise, reflect it about 
a vertical line and divide the weight of each right lozenge of the resulting region by q x+a to get 
the hexagon Hex(m,a,y) weighted by wt 2 - Then we get (4.14) from Corollary 3.2. The case 
t = 0, can be treated similarly to the case z = 0, based on Figure 4.2(d). 

The induction step is completely analogous to that of the proof of Proposition 4.1. We also 
apply Kuo’s Theorem 2.1, based on Figure 4.3. After removing lozenges forced by the shaded 
unit triangles, we get back new B-type regions weighted by iut 3 . Figure 4.3 tells us that the 
product of M 3 -generating functions of two regions on the top is equal to the product of the Ms- 
generating functions of two regions in the middle plus the product of M 3 -generating functions 
of two regions on the bottom. To precise, we get the following recurrence 

Ms (B mt a(x, y, z, t)) M 3 (B mt a(x, y-l,z,t-l)) = 

M 3 (B m ^ a (x, y-l,z, t)) M 3 (B m ^ a (x, y,z,t- 1)) 

+ q m+a+x+ y +z M 3 (B mt a(x, y — 1, 2 + 1, t — 1)) M 3 (B m ^ a (x, y,z- 1, t )), 

(4.16) 

and the proposition follows from the induction hypothesis. □ 

3 We still follow the process in Figures 4.2 and 4.3, however, the reader should be aware that the weight 
assignment here is different from that in the proof of Proposition 4.1. 
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x+m 



x+m y 


y+m 



Figure 5.1: The four base cases: (a) x = 0, (b) y = 0, (c) z = 0 (b), and t = 0. 


5 Proof of Theorem 1.2 

By Proposition 3.1, we only need to show that 


M 2 Q 


x y z t 


x y z t 


H q (m + a 


z + f) 


H q (m + a + b + c + x + y + t) H q (m + a + b + c + x + y + z) 
H q (m + a + b + c + x + t) H q {m + a + b + c + x + y) H q {m + a 


H q (m +1 


z + t) K q (m + a + b + c + x) H q (m + a + b + c + y) 


H 9 (x) H,(y) H q (z) H 9 (t) H,(m ) 3 H 9 (a ) 2 11,(6) H,(c) H q (m + a + b + c) 


H q {x + t) H q (y + z) 


H q (m + a ) 2 H q (m + b) H q (m + c) 


z + t) H q (m + a + c + x) H q (m + a + b + y) 
H q {m + b + y + z) H,(m + c + x + t) 

_H,(c + x + t) H,(6 + y + z) _ 

H,(a + c + x) H,(a + 6 + y) H,(6 + c + z + t) 


(5.1) 


We prove (5.1) by induction on y + z + t. The base cases here are the situations when one 
of the four parameters x, y, z, t equals to 0 . 
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We assume that our region Q := Q 


is weighted by wt 2 


x y z t 
K m a b Cj 

If x = 0, by applying Region-splitting Lemma 2.2, we split the into two parts as in Figure 
5.1(a): the shaded hexagon Hex(a, c,m) (weighted by wtz) and Q — Hex(a,c,m). We remove 
forced vertical lozenges (which have the weight 1 ) from the latter region, then rotate the resulting 
region 60° counter-clockwise, and reflect it about a vertical line. This way we get the magnet 
bar region B b>m (a,t + c,z,y) weighted by u>f 3 . Thus, we have 


“ 2 ( C (m l b 0) = “ 2 “3 l l ‘) 

= M 2 [Hex(a, c, m)) M 3 (B bim (a, t + c,z, y )), (5.2) 

and (5.1) follows from Corollary 3.2 and Proposition 4.2. The case t = 0 can be treated similarly 
to the case x = 0, based on Figure 5.1(d). The only difference is that our forced right lozenges 
have weight product equal to q x ( a + m )( z + b + c )+ x { + t + )_ Thus, we get 


M 2 Q 


x y z 0 
m a b c 

= M 2 (Hex(a, x + c, m)) M 2 (Q ( 1 ^ ) — Hex(a, x + c, m) 


m a b c 

= M 2 ( Hex(a , x + c, m)) q ^ + rn ){z+b+ c) + ^^ +1 ) m 3 {B b>m (a, c, 0 , y)) 

Again, (5.1) follows from the Corollary 3.2 and Proposition 4.2. 

If y = 0, again, by Region-splitting Lemma 2.2, we get 


(5.3) 


M 2 Q 


x 0 z t 
m a b c 


= M 2 ( Hex(m , b, a)) M 2 ( Q ( ^ * ) — Hex(m , b , a) 


(5.4) 


(see Figure 5.1(b)). We also remove forced lozenges (having weight 1) from the second region 
on the right-hand side to get a region R'. Next, we rotate R' 60° clockwise and reflect it about 
a vertical line to get the magnet bar B c , m (b + z, a, x, t ) weighted by wt\. Thus, we have 


M 2 Q 


0 z t 


m a 


= M 2 ( Hex(m , b , a)) Mi (B c>m (b + z, a, x, t )), 


(5.5) 


and (5.1) follows from Corollary 3.2 and Propositions 3.1 and 4.1. The case z = 0 can be 
obtained in the same way, based on Figure 5.1(c). 

For induction step, we assume that x,y,z,t are positive and that (5.1) holds for any Q-type 
regions in which the sum of the y-, z- and f-parameters strictly less than y + z + t. 

f x y z 

We now apply Kuo condensation to the dual graph G of the region Q I ' ) weighted 

by wt 2 - The four vertices u, v, w, s correspond to the four shaded unit triangles in Figure 5.2(b). 
By collecting the weights of the forced lozenges shown in Figures 5.2(b)—(f), we get respectively 

M(G-{«, v, w, a}) = q ( z+ T 1 )+(*+y+ m -2)(z+t+m+a+b+c) M 2 (q V ~ 1 £ t ^ ) , ( 5 . 6 ) 
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x+y+m 


x+y+m 



(a) 



x+y+m 



(c) 



(d) 
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x y — 1 z t 
m a be 


M {G - {u,v}) = q( Z+ ™ +1 ) M 2 (q 
M(G - {w, a}) = q ^+y+m-2){z+t+m+a+b+c) ^ j^g ^ 

M(G - {u, a}) = g( Z+ 2 +1 ) M 2 (q ( 


x y z t — 1 

m a b c 


x y — 1 z + 1 t — 1 
m a b c 


and 


M (G - {v, ttf}) = q ^+y+m-lKz+t+m+a+b+c) ^ ^g ^ 

Substituting (5.6)-(5.10) into equation (2.1) in Kuo’s Theorem 2.1, we get 


x y z — 1 t 
m a b c 


(5.7) 

(5.8) 

(5.9) 

(5.10) 


m 2 (q( x y i MWgf* y ~ l i t - 1 

' \m a b c J \ \m a b c 


"a 1 1 0W< 11 A 


+ z+t+m+a+6+c Ma | g / a; y 1 2 + 1 t 2/ z 1 f 

1 a b c JJ \ \m a b c 


(5.11) 


/x yy z t\ 

Finally, if we denote by T ( ’ J the expression on the right-hand side of (5.1), 


we 


only need to show that T satisfies also the recurrence (5.11). Equivalently, we need to verify 
that 

vi, ( x y - 1 2 A n,( x y z l ~ l 

\m a be) \m a b c 


M/| x V z t) x y-l z t-1 

,m a b c) \m a b c 


q 


z-\-t-\-m-\-a-\-b-\-c 


^ | x y z-1 A^/x y-l 2 + 1 t-1 
m a b c) \m a b c 




fx y 

z 

A 

+ ( x y^ 1 z 

\m a 

b 

c) 

\m a b c ) 


= 1. (5.12) 


Let <L ( X y Z A ; = q \ m a b c ) f x y ~ ). By the definition of the func- 

\m a b cj \m a be) 

tion g, we get 

(x y z A, ( x y z i —1\ (x y z A, lx y — 1 z t — 1\ 

S \m a b c) +S [rn a b c J ~ S [m a b c) + g \m a b c ) ( ' 5 ' 13 ^ 

and 


x y z t 


x y 2—1 A lx y 2—1 t 

S a b c) \ m a b c 


= (m + x + y- z-1)+ g( X V ~ A + g f 

\m a b c) y 


x y — 1 2 t — 1 

m a b c 


■ (5.14) 
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Therefore, ( 5 . 12 ) is equivalent to 


x z/ — 1 z t 


7 $, , 

,m a be \m a b c 


x y z t — 1 


$ 


x y z t 
m a b c 


$ 


x y — 1 z t — 1 


m a 


m-\-a-\-b-\-c-\-x-\-y-\-t —1 


$ l s y z-l A $ / x y -1 2 + 1 t-1 

m a b cl l to a b c 




x y z t 


$ 


X y — 1 2 t — 1 


= 1. 


(5.15) 


Let us simplify the first term on the left-hand side of ( 5 . 15 ). We notice that the two $- 
functions in the numerator and the denominator of the first fraction in the first term are different 
only at their y-parameters. Canceling out all terms, which have no y-parameter, and using the 
trivial fact H q (n + 1)/H q (n) = [n] 9 !, we get 


$ 


x y — 1 2 t 
to a be 


$ 


x y z t 


[y + 2 - l] g ![a + b + y - l] g ![m + b + y + z - l] g ! 
[y - 1 ]<?![& + V + z - 1 } q \[m + a + b + y- 1] 9 ! 


to a b cj 

[■m + a + b + c + y — l] g ![m + a + b + c + x + y + t — 1] 9 ![to +a+b+c+x+y+z— l] g ! 

[m + a + b + c + x + y — l] g ![m +a+6+c+y+2— 1 ] 9 ![to + a + b + c + x + y + z + t — \] q \ 

(5.16) 

Doing similarly for the second fraction of the first term, we obtain 
' x y z t — l\ 

m a b c ) _ [V - !]?![& + V + z- l] q \[m +a+b+y- 1],! 


$ 


$ 


x y- 1 2 t- l\ [y + z - l] g ![a + b + y - l] q l[m + b + y + 2 - l] g ! 


m a b c 

{m + a + b + c + x + y — l] ? ![ro + a + b + c + y + z — l] g ![ro + a + b + c + x + y + z + t — 2 ] q \ 


[m + a+b + c + y— l] g ![ro + a + b + c + x + y + t — 2 ] q \[m + a + b + c + x + y + z— l] q \ 


This implies that the first term on the left-hand side of ( 5 . 15 ) can be simplified as 

[m + a + b + c + x + y + t — l] q 


+ i x y - 1 z $( x y z t ~ l 

, to a be) \ to a b c 


^ , x y z A /a; V~ 1 z A [TO + a + & + c + x + y + 2 + t-l] g 

,to a b cl \ to a b c 


(5.17) 


(5.18) 


We simplify the second term on the left-hand side of ( 5 . 15 ) in the same way (the numerator 
and the denominator in each fraction are now different at their 2-parameters). We get 


$ |X y 2-1 A $ /x y — 1 2 + 1 t- 1 

.too b cl l to a b c 


[4 


$ 


x y z t 


$ 


x y— 1 2 t — l\ [m + a + b + c + x + y + z + t—l] q 


(5.19) 


.to a b cl V to a b c 

By ( 5 . 18 ) and ( 5 . 19 ), the equality ( 5 . 15 ) becomes the following identity 

[m + a + b + c + x + y + t- l\ q + q ™+a+b+c+x+ y +t- 1^ 

[m + a + b + c + x + y + z + t— l] q 
which follows directly from the definition of the q- integers. This completes our proof. 


= 1 , 


(5.20) 
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